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Section I 

INTRODUCTION 


In launching the Space Shuttle, it is expected that two recoverable Solid 
Rocket Boosters (SRB) wi’l be used. These boosters will detach from the main 
vehicle and drop to the ocean, from which they are to be recovered. In some 
cases, a component, or possibly a whole booster, will not be recoverable and 
will have to be replaced. 

It is important to know the probable number of booster components needed 
to successfully accomplish the Space Shuttle Mission of 440 flights, so that 
advance planning on the optimum number of SRB components to be built can be 
carried out. Several situations are analyzed, namely, these cases where: 

(1) The booster is used continuously until not recoverable; 

(2) The booster is used for 20 launches unless not recoverable; 

(3) The booster is made up of identical. Independent, replaceable compo- 
nents . 

Two independent techniques have been used in this analysis. The first 
method, to be called the State Probability Method, consists of defining an 
appropriate state space for the outcome of the random trials. Then the proba- 
bility of reaching each state is evaluated. The second method is called the 
^Monte Carlo Technique, or the Model Simulation Method. Detailed discussion of 
the Monte Carlo Method can be found in references 1 and 2. The Monte Carlo 
Method is simpler to formulate than the State Probability Method, but requires 
longer computer running time and can have statistical uncertainty due to sampling 
size. Both methods are compared to check the results. 

Finally, a simple analysis of expected cost is carried out. This includes 
the case where the boosters can be built initially at some initial cost and the 
remainder built as needed. The resulting optimum number to build initially and 
the minimum expected cost are evaluated. 

1. Graves, M. E. and Perlmutter, M. "Statistical Analysis of Flight Times 
for Sjpaos Shuttle Ferry Flights ", MSA CP-129021, Jan 1974. 

2. Shreider, Y. U. A., m . The Monte Carlo Method ", Tergamon Press, 1966. 



Section II 


STATE PROBABILITY AND MONTE CARLO METHODS 

The State Probability Method consists of defining the statistical problem 
to be solved in terms of a state space of pertinent possible outcomes of the 
random trials. The probability of obtaining each of these outcomes is then 
calculated. The method will be illustrated in the solution to the present 
problem. 

The Monte Carlo Method simulates the problem, using a random number gener- 
ator on a computer. A sufficient number of trials are simulated to obtain the 
statistical behavior of the system. 

The probability density of the number of one-component boosters needed for 
a mission of 440 launches is calculated in subsections 2.1 and 2.2. This is 
found for a given attrition rate of losing the booster in a launch. The boos- 
ter is used continuously until nonusable. 

Subsection 2.1 gives the State Probability Method of solution, while sub- 
section 2.2 gives the Monte Carlo solution. 

Subsection 2.3 and 2.4 consider the case where the booster can only be 
used for a maximum of 20 launches. Subsection 2.3 considers the State Proba— 
bilty solution while subsection 2.4 uses the Monte Carlo solution. 

Subsections 2.5 and 2.6 extend the results to the multicomponent rocket 
booster. In this case, the booster is made of Identical replaceable components. 
Each component has the same attrition rate, q. The probability distribution 
of the number of components needed for 440 missions is calculated. In subsec- 
tion 2.5 the State Probability Method is used, while in subsection 2.6 the 
Monte Carlo analysis Is used. 

2.1 STATE PROBABILITY METHOD - CONTINUOUS USB, ONE-COMPONENT BOOSTER 

To evaluate the probability density of the number of one-component boosters 
needed to achieve a mission of 440 launches for a given attrition rate, q, per 


2-1 


launch where the booster is used continually until lost, then replaced, the 
following steps are taken* First, a state is defined as shown in Figure 

2-1. The term S,.^ refers to the state of k launches with r booster replace- 
ments. The state probability for the booster replacement in 1 launch is given 
by 


P^l.k* * Pk 


(2-1) 


where is the probability of replacing one booster in k launches* 



* 




s h DEFINES NUMBER OF BOOSTERS REPLACED, b, DURING A TOTAL NUMBER 
’OF LAUNCHES, f 

P f DEFINES THE PROBABILITY OF FAILURE OF ONE BOOSTER IN f LAUNCHES 



Figure 2-1. STATE SPACE DEFINING NUMBER OF BOOSTERS REPLACED FOR TOTAL NUMBER 
OF LAUNCHES 



The attrition rate, or probability of booster loss per launch, is q, so 
letting p-l-q, then 


p 2 »qp 

k-1 

pk-qp* 


( 2 - 2 ) 


Therefore, 


I p (S. k ) * <1 [1 + p + p 2 + . . . ] « 3 — . i 
k-1 ’ 1_p 


(2-3) 


The probability of 0 booster replacements in the mission of 440 flights is 
given by 

439 


»io-I p(s x ,) - l-q 1 p^" 1 - p 439 

1,0 j-440 1,J j-1 


(2-4) 


This is the same result given by the binomial distribution for 0 replacements 
in 439 launches. 


P x Q - B (k - 439, r - 0) 


. v kl k-r r 

where, B(k,r) - r f ( k 1 r ) ~ l P q 


(2-5) 

( 2 - 6 ) 


This process is continued to the state S 2 which refers to replacing two 
boosters in k total launches. The probability of reaching state S 21t as can 
be seen from Figure 2-1 is given by 


P <s 2|J ) 


P < s l,j-i> P x + P (S l,j-2 ) p 2 + **• + P (S i,i ) P J-1 

2p J-l p l + 2p J-2 p 2 + - + 2p J±l P Jzl J ° dd „ 

2 2 1 

2p J-l P 1 + 2p j-2 p 2 + •** + P i P i • J * V *“ 

2 2 


7 ) 
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This can be seen to be equal to 


p( s 2f j> a (j-i) q 2 p^~ 2 


(2-8) 


This can be checked as before 


l P(S, ,) - l (3-D q 2 pJ- 2 - q 2 [1 + 2p + 3p 2 + ... 1 
j-2 ,J j-2 


(2-9) 


(1-p) 2 


The probability of using one or fewer booster replacements in the mission is 
then given by 


• 439 439 

* i » l P(S 2 ,) - 1 " l P(S, .) - 1 - l (3-1) q 2 pJ" 2 • 
* 3-440 ’ J 3-2 j-2 


(2-10) 


It can be seen that (ref. 3, Section III) 


. _ 

I (3“1) P^~ 2 - 1 + 2p^ + 3p 2 + , . . + 438p 
3-2 

. 1 + p (ll£i 37 _> . U +_437J p' 


( 2 - 11 ) 


Then, 


1 - *38 (-1 + * 38 n 

l 2 q 2 « 


P - P 438 (q + P) + 438q p* 38 - p 439 + 439q p 43 ® 
* »1 


(2-12) 


This is equivalent to the cumulative binomial distribution with 1 or leas 
replacement* in 439 launches (see eq. 2-6) 


p i l * " *39, r ■ 0) + B (k ■ 439, r - 1) 


(2-13) 


8 . Jolley, L. B. W. t "Stemation of Soviet", Dover Publioatine 1961. 
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This process can be continued to give 


P 1 2 * B(k ' 439 » r - 0) + B(k - 439, r - 1) + B (k - 439, r - 2), 
etc. ’ (2-14) 

Besides the cumulative distribution P. that is, the probability of 

i , r 

using r replacements or less in the mission, p. can be written, the distri- 

i t r 

bution which gives the probability of using exactly r replacement- for the 
mission. Thus, 


P 1 0 " B r * 

p,^ - B (k - 439, r - 1) (2-15) 

p ■ B (k » 439, r • 2) 

1,2 

The probability distribution Pi >n shown in Figure 2-2 gives the proba- 
bility of using n one-component boosters in a mission of 440 launches. 
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The cumulative distribution, Pi >n which gives the probability of using n 
or less one-component boosters in 440 launches, is shown in Figure 2-3. As 
has been shown previously, the probability density in Figure ’ - ,'s the 

binomial distribution for 439 launches or trials where th. probabi. ty of 
failure or attrition is q and the number of failures, or boosters replaced, is 
r The mean and variance of the binomial is given by 



U - 439q and a ^ ■ 439q (1-q). 


rrrTT - T , | ,in 

COMPONENTS PER 
EQUIVALENT BOOSTER 

1 - COMPONENT BOOSTER I 

1 - COMPONENT GAUSSIAN APPRO X. ? 

• 1 - COMPONENT MONTE CARLO 1 

O 2 - COMPONENT BOOSTER 2 

+ Q - COMPONENT BOOSTER 8 


(2-16) 
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Figure 2-3. CUMULATIVE PROBABILITY OF USING n OR LESS COMPONENTS IN 1, 2 
OR 8 COMPONENT BOOSTERS IN 440 LAUNCH MISSIONS; BOOSTERS RE- 
PLACED AT FAILURE 
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For a large number of launches, the binomial distribution can be approxi- 
mated by a Gaussian distribution with the same mean and variance as the 
binomial. The 95 percent probability of using r replacements of the booster 
in 439 launch missions is given by the Gaussian distribution as 


1.654 - 


(2-17) 


The 50 percent (Pj^ n » .5) and 95 percent (Pi #n ■ «9 e ) cumulative probability 
for number of single-component boosters needed for a mission of 440 launches 
is plotted in Figure 2-4. 




2.2 MONTE CARLO METHOD - CONTINUOUS USAGE, ONE-COMPONENT BOOSTER 

For the same conditions discussed in the previous subsection a Monte Carlo 
solution can also be obtained. A set of pseudo random numbers which are 
uniformly distributed between 0 -nd 1 is generated. Sequences of these numbers 
are examined and recovery of SRB is assumed if the values are less than 1-q 
and non-recovery of SRB is assumed for any number exceeding 1-q. The probability 
of 440 successful missions is calculated for a given number of starting SRB. 

These results are plotted in Figure 2-3 and are in good agreement with the 
results of thi State Probability Method. 


2.3 STATE PROBABILITY METHOD - 20 USES ONE-COMPONENT BOOSTER 

A calculation similar to that of the continuous usage case can be carried 
out if the booster is replaced after 20 uses or failure. For this case, Instead 
of (eq* 2-2), the probability of replacing one booster in k launches is obtained. 


p k - q p k "\ 1 < 4 < 20 

p 20 - a>p 19 


(2-18) 


P k - 0 ; l > 20 
As before, (eq. 2-1) gives 


XS 1>k > - P k . (2-l») 

How it esn be seen thet 


l PC®. J - l P, - qll ♦ P ♦ P 2 ♦ ... ♦ P 18 ] + P 19 
k-1 A »* j-1 J (2-20) 

2 17 18 

■ q(l p + p + p A/ ) ■* p AO ■ 1 


Then es before. 

•* J “ 

- I P(s. k ) . (2-21) 

k*l 


1.0 


I p<s ltk ) - 1 


k-440 
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This process is continued as before to yield 


1J7 

f i.l ‘ 1 • & P<S M ) 1 


( 2 - 22 ) 


where p(S 0 is given by equation 2-7. This procedure is then continued to 
obtain P^, ? 1 3 ... . 

The probability distribution for the number of single-component boosters 

for 20 uses in a 440 launch mission is shown in Figure 2-2, and the cumulative 

distribution which gives the probability that n or less booster components will 

be used for the mission is shown in Figure 2-5. The result can be seen to be 

close to Gaussian, since the cumulative distribution is close to a straight line 

on the probability graph in Figure 2-5. The 50 percent (p, =0.5) and 95 percent 

i,n 

(p, =0.95) cumulative probability for the number of boosters needed are also 

1 , n 

plotted in Figure 2-4, along with the case for continuous usage of the boosters. 
It can be seen that for high attrition rate, q = 0.2. The results for these 
two cases converge. 

2.4 MONTE CARLO METHOD - 20 USES, ONE-COMPONENT BOOSTER 

The Monte Carlo procedure in the present case of a maximum 20 uses of the 
booster is similar to the Monte Carlo solution for the continuous use case, 
except in the simulation the booster is replaced after 20 consecutive launches. 
Again the probability of the number of boosters needed to achieve 440 launches 
is calculated and plotted in Figures 2-2 and 2-5. 

2.5 STATE PROBABILITY METHOD - MULTI-COMPONENT BOOSTER 

The SRB can be assumed to be built of replaceable components, each 

component having an independent attrition rate. Interest is in the probable 

total number of components needed for a mission of 440 flights if each component 

has Cite same attrition rate, q. Again, a state space c. is drawn as shown in 

1 > f 

Figure 2-6. The state c. represents 1 Independent components having a total 

i 1 1 

of r replacements during a mission of 440 launches. 

Considering a one-component rocket: 


p(ci >r ) ■ P 1>r 5 r = 0, .... 439 


(2-?3) 


2-9 



NUMBER OP EQUIVALENT BOOSTERS USED 















+ 










_e 



TTT" l'~T 7 ITTTT 

COMPONENTS PER 
EQUIVALENT BOOSTER 

I - COMPONENT BOOSTER 1 

] - COMPONENT GAUSSIAN APPROX. 1 

• 1 - COMPONENT MONTE CARLO 1 

O ? - COMPONENT BOOSTER Z 

4 2 - COMPONENT MONTE CARLO ? 

+ 8 - COMPONENT BOOSTER 8 

* 8 - COMPONENT MONTE CARLO 8 

i — '1 — i — r i—i — M — I — h 


q = ATTRITION RATE 
FOR COMPONENT PEP LAUN 



MT+ . 

i 1 

u. -f- 

i 

t — r 

i . 


q = . 15 - 


+ + 


43 *- 

I * ! 

-4— -t * 


I* ^ 


q * 0.04 



Figure 2-5. CUMULATIVE PROBABILITY OF USING 1 OR LESS COMPONENTS IN 1, 2 OR 
8 COMPONENT BOOSTER IN 440 LAUNCHES: BOOSTER REPLACED AFTER 20 
USES OR FAILURE 





Figure 2-6. STATE SPACE DEFINING THE TOTAL NUMBER OF COMPONENTS REPLACED IN 
THE MISSION FOR A ROCKET BOOSTER BUILT WITH A GIVEN NUMBER OF 
REPLACEABLE COMPONENTS; C c ,r DEFINES NUMBER OF COMPONENTS REPLACED, 
r DURING A SUCCESSFUL MISSION USING ROCKET BOOSTER CONTAINING c 
REPLACEABLE COMPONENTS; P c , r IS THE PROBABILITY OF STATE C c ,r 
OCCURRING 



where p is the probability of r replacements of the one-component during a 
I , r 

successful mission. This result was obtained previously in equation 2-15, etc., 
and is shown in Figure 2-2. 


Obviously, 


r»0 


■l.r 


(2-24) 


As seen in Figure 2-6, state which represents the state of a two- 

component booster with a total of r component replacements in the mission can 
be found as follows. Consider the state c^ q, that is a two-component booster 
with zero component replacements in the mission. The probability of obtaining 
this state is given by 


p<c 2,o ) ■ p 2,o ■ ric i.a> p i,o ' p i,o p i,o • 

This process can be continued to give 

' P 2,l ' P<C 1.0 ) p l,l + p(c l,l ) p l,0 " 2p l,l P l,0 

or for an odd number replacements, r odd, 


(2-25) 


(2-26) 


p 2,r " ’ 2p l,r P l,0 + 2p l,r-l ,p l,l + — 

i, 2 L * 2 


(2-27) 


while for r even, 


P 


2,r 




• • • 



(2-28) 


This can be extended to the state c. , which represents a 4-component booster 

*>»r 

with a total of r components replaced in the mission as shown in Figure 2-6. 
Continuing as in the previous case, the probability of achieving state c^ Q , 
that is, a 4-component booster with 0 component replacements in the mission, 
can be written as 
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P 2,0 P 2,0 


(2-29) 


p 4,0 ' p < c 4,0> ‘ ' ,(c 2,0 ) p 2,0 


Then continuing as before, 


P 4,l = P<C 4,1 ) " p(c 2,0 ) P 2,l + p(c 2,l ) P 2 ,0 “ 2p 2,0 P 2,l 


P 4,2 " 2p 2 ,2 P 2 ,0 + P 2, 1 P 2,l 


(2-30) 

(2-31) 


This can be extended for r odd to 


p 4,r ’ p(c 4,r> ’ 2p 2.r p 2,0 + 2p 2,r-l p 2,l + "• + 2p rH ", r4 

1 ’ 2 2 


(2-32) 


and for r even. 


p 4,r- 2p 2,r p 2,0 + - +p r» r 

* 1 2 ^ »2 


(2-33) 


This procedure can be applied to 8 similar components in a similar manner, giving 


P S,0 ‘ p<c 8,0> ■ P(C 4,0> p 4,0 ■ p 4,0 p 4,0 


(2-14) 


p 8,l ■ p(c 8,l ) ■ P(C 4,1 ) P 4,0 + P(C 4,0> p 4,l ’ 2p 4,l p 4,0 (2 - 35) 


Again, for a total of r replacements of an 8 independent component booster for 
r odd. 


p e.r ■ 2p 4,r p 4,0 + 2p 4,t-l p 4,l + — ' 2p , £±1 P lJz l < 2 - M > 

while for r even. 


4 4 - 

2 2 


2Px 


+ ... + p r P . 


r 8,r 4,r **4,0 

These reaulta are shown in Figures 2-3 and 2-5 as cumulative distributions. 


(2-37) 


A — A — 

4, 2 4, 2 


; - r ' Jo 


(2-38) 
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2.6 MONTE CARLO METHOD - MULTI-COMPONENT BOOSTER 


The simulation for the multi-component SRB is similar to the Monte Carlo 
analyses for the single-component SRB. Two sets of pseudo random numbers are 
used for the 2-component booster, while 8 sets of pseudo random numbers are used 
for the 8-component SRB. The same method as used in the single -component SRB 
is then used. The results of these calculations are shown in Figure 2-5. The 
probability of having a specified number of booster components remaining after 
440 missions, given the number of booster components initially available, is also 
shown in Figure 2-7. 


The confidence range of the present Monte Carlo calculation can be found 
following the procedure shown in reference 1. Assume Mg is the number of 440 
launch missions using n or less boosters in N Monte Carlo trials. Then che 
probability of this occurring p £ R , is given as 


M 

f Z p c,n 


(2-39) 


M 

The 99 percent confidence ranges, (see ref. 1) on is given by 


M 

s _ 

N ” P c,n 


6 < 2.57 /p„ (l-p r ) /N 
— c,n c,n 


(2-40) 


If, considering p ■ 0.98 one may obtain for N-350 trails 

C|ll 

6 < 2.57 [°-y)] - 0.019 . (2-41) 

This result gives an error range in p of +0.019. This would correspond, 

c f n mm 

using the two-component State Probability results, to an equivalent error range 
of +6 equivalent boosters or +12 components. The Monte Carlo results fall 
well within this range, as can be seen in Figure 2-5. 
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Figure 2-7. PROBABILITY OF k OR LESS BOOSTER COMPONENTS REMAINING AFTER 440 

MISSIONS GIVEN n INITIAL BOOSTER COMPONENTS WITH THE ATTRITION RATE 
q OF 0.04 (LIMIT 20 REUSES OF EACH COMPONENT); TWO COMPONENT BOOSTER 





Section III 

APPLICATION OF RESULTS TO SHUTTLE SRB COST ANALYSES 


The present results can indicate expected costs of the SRB and aid in the 
decision-making process to achieve minimum costs. 

If the cost of building a booster is assumed as c, then the expected cost 
of building the boosters for the 440 launch mission is given by 

00 

E.C. » C y np - c u (3-1) 

£ r c f n n 

n=l 

This shows that the expected cost of the booster for the mission is equal to the 
average of the distribution of number of boosters used in the mission times the 
cost per booster. Figure 2-4 indicates that the average number of boosters for 
the mission, which is approximately equal to n * 0.5, and also therefore 
the cost, will rise with the attrition rate and is higher for the 20 reuses 
case than for the continuous usage case. If we assume the multi-component 
booster cost is equal to the cost of a one-component booster, then, since the 
mean for the multicomponent booster is close to the value for the mean of a one- 
component booster (see Figures 2-3, 2-5), the cost for the multi -component 
boosters should be close to the cost for the single -component booster. How- 
ever, the components are likely to have a significantly lower attrition rate 
which would give significantly lower cost for the multi -component booster 
compared to the single -component booster. 

The decision process can be extended to the case trhere an initial group of 
boosters can be built at one cost and the remainder built as needed at a 
dilferent cost. An objective is to find the optimum number to build initially 
so as to obtain a minimum cost. This can be done as follows. Assume an 
initial cost of building booster components of C^. Then if additional booster 
components are needed to complete the mission, they can be built at cost c f . 

A decision table can be constructed as shown in Table 3-1. 
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Table 3-1. DECISION TABLE FOR BOOSTER COMPONENTS 
COMPONENTS BUILT INITIALLY 


oc 

© 


Q 

LU 

Q 


LO 

► 

ZT 

LU 

Z 

o 

Q_ 

z: 

o 

o 


to 

s 

CO 


O 

OC 


NUMBER OF COMPONENTS BUILT INITIALLY 
b* 1 2 


"l.l 

c <»l,l>- c 1 

EC<M,i). c iH u ,)P c ., 

"1.2 

C l"l.2»- 2C 1 

ec(» 1i2 K(», i2 )p C i , 

"1.3 

C(N l,3 )=3C i 

EC(N 1>3 )= c (N lt 3)p c> i 

1 

N 2,l 

c(N 2t1 K. + c f 

EC(N 2f iK(N 2>1 )p c>2 

N 2,2 

C(N2 i2 )=2 Ci 

EC(N 2f2 )= c (N 2> 2)p Cf2 

N 2,3 

c (N 2f3 )= 3C 1 

EC(N 2 , 3 )- c ( N 2,3)Pc.2 

1 

N 3,l 

C < N 3,l )=C 1 +2C f 

EC(N3 fl K(H 3i1 )p Ct 3 

N 3,2 

C{N 3,2 ) ‘ 2C i +C f 

EC(N 3f2 ).C(N3 >2 )p c> 3 

N 3,3 

C(N 3 , 3 > = 3C 1 

EC < N 3.3 )<(N 3,3 )p c,3 

1 

• 

i ♦ 

• 

t 

• 

1 

«»• j«V.»c 

EC(2). 

EC(3) ' n I, C(N o.3 )p c.n 

1 


The possible decision on the nuaber of coaponents to build Initially is 
along the top of the table and the nuaber needed for the aission along the aide. 
The state N . ia defined, where b refers to the nuaber of components built ini- 

Q p 0 

tially and n refers to the number needed* The cost of each state is given by 

C(N . ) • bC. + (n-b)C. for n > b (3“ 2 ) 

n,o i i 

■ bC^ for n < b 
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This can be normalized, 


C(N n h> 

n,D 


C(N n h> 
n » b . 

C i 


b + (n-b)(C £ /C £ ) for n > b 
b for n < b 


(3-3) 


The probability of each state occurring can be written as 


p(N 



P 


c,n 


(3-4) 


where the subscript c refers to the number of components per booster and n 

refers to the number of components used per mission. Then the expected 

cost of each state is given by EC(N , ) • c(N . )p . The expected cost of 

n,D n,D c,n 

each decision can then be found to be given by 

00 

EC(b) - I c(H ), . (3-5) 

n-1 * 

The optimum decision of the number of booster components to be built initially, 

b , is the value that causes the expected cost to be a minimum EC so that 
n m 

EC(b ) - EC . (3-6) 

n o 


These results have been evaluated which are shown in Figures 3-1 and 3-2 as a 
function of the cost ratio C £ /C £ . 

The results in Figure 3-1 are plotted in percent cumulative probability of 
the number of booster components used in the mission. The results indicate that 
if the ratio of final cost per booster to initial cost C £ /C £ is less then some 
minimum value, the optimum procedure is to build one booster initially, then 
build the rest ss needed. For values, of C £ /C i greater then this minimum, 
build initially most of the components needed for the mission so that one is 
80 percent or more confident that sufficient boosters are buiJt to complete the 
mission, even though some may be left over. 


In Figure 3-1 the normalized cost of the optimum decision is shown. Itach 
more detailed calculations than those shown here are necessary for a complete 
cost analysis. ,..<1 


CUMULATIVE PROBABILITY OF NUMBER OF COMPONENTS USED IN MISSION 












Section IV 
CONCLUSIONS 


The Space Shuttle wi 1 recoverable solid rocket boosters. However, at 
tines due to damage or wv^r the rocket booster will not be recoverable and must 
be replaced. It is oi ’ -.^ rr.ance to know the probable number of solid rocket 
boosters chat will ; . a oe replaced in the nominal mission of 440 Space 
Shuttle launches. information is necessary to the planning of the mission 

so that the solid rocket booster replacements will be available as needed at 
a minimum cost. 

The present analysis calculated the probability distribution of solid 
rocket boosters needed for a mission of 440 launches for a given attrition rate 
or probable loss of the solid rocket booster per launch. Two methods of cal*- 
culatlng the probability distributions were used. Ihese are the State Probability 
Method and the Monte Carlo Method. The results from both methods agreed. It 
was found that the Monte Carlo Method was easier to formulate, while the State 
Probability Method required less computing time and was more accurate. 

The results of the calculation showed that the probability distribution 
of booster replacements for the 440-launch Shuttle mission was close to 
gaussian with a mean value that rose linearly with the attrition rate. Thus, 
as can be seen in Figure 2-4, for an attrition rata of 0.04 the mean or expected 
value of boosters needed for the 440-launch mission, which in this case is 
equivalent to the 50 percent cumulative p obabl • '.ty, is 18. This number goes 
up to 44 boosters for an attrition rate of 0.1 and Increases further to 88 
boosters for an attrition rate of 0.2. 

Also shown in Figure 2-4 is the number of boosters for the 95 percent 
cumulative probability value. The probability is 0.95 that this nuaber of 
boosters or leas will be needed to successfully complete the mission of 440 
launches. As shown in Figure 2-4, at an attrition rata of 0.04 the 95 percent 
cumulative probability nuaber is 25 boosters. This nuaber rises to 54 boosters 
for an attrition rate of 0.1, and to 103 boosters at an attrition rate of 0.2. 
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The effect of wear on the solid rocket boosteis was included in the 
analysis by replacing the booster alter 20 launches. This caused the expected 
number of rocket boosters needed to complete the mission to increase. As can 
be seen from Figure 2-4, the mean or expected values of solid rocket boosters 
needed for a 440-launch mission changed from 18 to 33 for an attrition rate of 
q ■ 0.04 if the rocket booster was replaced after 20 launches instead of being 
used continuously. At an attrition rate of 0.1 whe number of boosters needed 
went frotr 44 to 49 in going from continuous use to 20 uses, while at an at' it ion 
rate of 0.2 the number of. boosters needed was the same for the case of 20 uses 
as for the continuous use case. 

The boosters can be built up in replaceable components. If damaged, these 
components can be replaced without replacing the entire booster. An analysis 
was carried out to calculate the probability istribution of the number of 
replaceable components needed for a 440-launch mission for a given component 
attrition rate. This result is shown in Figure 2-3 for the continuous booster 
use case. As can be seen from the fif re the expected number of components 
needed for the mission is given in terms of equivalent boosters. An equivalent 
booster is equal to the number of components in a booster, so that 10 equivalent 
boosters for an eight component booster consists of eighty components. Thus for 
a component attrition rate of 0.04 it can be seen that the expected number of 
equivalent boosters for a booster msde up of eight replaceable components is 
nearly the same as thst for a booster with one component, end is 18. This would 
make the expected number of componente for the eight component booster equal 
to 144. 

For the case where the components are replaced efter 20 uses Figure 2-5 
is used. If the results are examined for a component attrition rate of 0.1 
we see that the expected number of equivalent boosters for the eight component 
booster is larger (52) than the expected number for the one component booster 
(50). This would make the expected number of components needed for the 440- 
launch mission with an eight component booster with an attrition rate per 
component per launch of 0.1 equal to 416. 
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A simple cost analysis was carried out, based on the probability distribu- 
tion obtained. The results indicate that the expected cost should be proportional 
to the expected number of boosters needed. Thus the cost would rise as the 
attrition rate rises and should be higher for the 20 reuse case than for the 
continuous use case; see Figure 2-4. Assuming the cost of the multicomponent 
equivalent booster is close to the cost of a single component booster, and since 
the expected values are close for both cases for the same attrition rate, the 
costs should be similar for the two cases. The attrition rate should be lower, 
however, for the multicomponent booster. This is because if one component in 
an eight component booster is damaged only one component needs to be replaced 
whereas for a one component rocket the entire rocket must be replaced. This 
would give an attrition rate for the eight component booster 1/8 that for the 
one component booster. This would result in an eight fold cost cut for the 
multicomponent booster, compared to the one component booster. 

The cost analysis was extended to the case where an initial group of 
boosters could be built at some initial cost and the remainder built as needed 
at a final cost. The objective was to find the optimum number to build initially 
as a function of the ratio of initial to final cost so as to minimize the overall 
cost. In Figure 3-1 the results show that if the ratio of final cost to initial 
cost is below a certain value the boosters are built as needed. However, if 
they are above this value, then enough boosters should be built initially so 
that one is fairly certain of having sufficient boosters to complete the 
440-launch mission. Thus from Figure 3-1 for a one component booster in 
continuous use with an attrition rate of 0.04, if the ratio of final to 
initial cost is below 1.54, build the boosters as needed; if this ratio is 
higher than 1.54, then build enough boosters so that the probability is 0.94 
that there are sufficient boosters to complete the mission. 

The analysis can readily be extended to other values of attrition rate, 
number of launches per mission, and number of components per rocket. 
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